Existence and nonlinear stability of steady 
states of the Schrodinger-Poisson system 

(N ; 

Peter A. Markowich*, Gerhard Rein^, Gershon Wolansky* 
i— j ; * Mathematisches Institut der Universitat Wien 

Boltzmanngasse 9, 1090 Vienna, Austria 
' International Erwin Schrodinger Institute 
> Boltzmanngasse 9, 1090 Vienna, Austria 

CM ' * Technion, 32000 Haifa, Israel 

O ■ 

o : 

Abstract 

We consider the Schrodinger-Poisson system in the attractive 
Q^' (plasma physics) Coulomb case. Given a steady state from a cer- 

r~| ■ tain class we prove its nonlinear stability, using an appropriately de- 

, fined energy- Casimir functional as Lyapunov function. To obtain such 

steady states we start with a given Casimir functional and construct 
a new functional which is in some sense dual to the corresponding 
• i— i . energy-Casimir functional. This dual functional has a unique maxi- 

^ | mizer which is a steady state of the Schrodinger-Poisson system and 

lies in the stability class. The steady states are parametrized by the 
equation of state, giving the occupation probabilities of the quantum 
states as a strictly decreasing function of their energy levels. 
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1 Introduction 



A large ensemble of charged quantum particles interacting only by the electro- 
static field which they create collectively can be described by the Schrodinger- 
Poisson system: 

i^ = -A^ k +V^ k , (1.1) 
Ay— n, (1.2) 

oo 

n = Y,h\i>k\ 2 - (1.3) 
fc=i 

Here ipk = 4 ! k{t,x) is the wave function of the kth state, fee IN, Afc>0 denote 
the corresponding occupation probabilities normalized such that X)fcAfc = l, 
n = n(t,x) is the number density, and V = V(t,x) the self-consistent potential 
of the ensemble. In order to avoid continuous spectra we shall analyze this 
system on a bounded domain Q C IR 3 with sufficiently smooth boundary, and 
we supplement it with Dirichlet boundary conditions: 

fa(t,x)=0, V(t,x) = 0, t>0, xedQ, fee IN. (1.4) 

We could also consider the system on the whole space IR 3 and add to V an 
appropriate exterior potential V e . Initial data are given by a complete or- 
thonormal system (ijj k (-,0)) in L 2 (Q). We refer to [fll 0], E| |TB| for background 



information on the Schrodinger-Poisson system ( |1 . 1|) , (|1.2|) , (jL 

In terms of the density operator R(t) of the system, a time dependent, 
hermitian, positive trace-class operator acting on the Hilbert space L 2 (Q), 
the time evolution is given by the von-Neumann-Heisenberg equation 

3D 

Here the Hamiltonian is defined as Hy := — A+V(t,x) with potential V 
given as the solution of the Poisson equation ( |1.2| ) with Dirichlet bound- 
ary condition, and n(t,x) : = p(t,x,x) where p(t,x,y) is the kernel of the op- 
erator R(t). The Schrodinger-Poisson picture and the Heisenberg picture 
are equivalent: Let (<p k ) be a complete orthonormal sequence of eigenvec- 
tors of R(0) with eigenvalues (A&) and let (i/; k (t,-)) be the solution of the 
Schrodinger-Poisson system (|1.1|)-(|L~4|) with initial data tpk(0) — <Pk- Then 
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p(t,x,y) =J2k^k"ipk(t,x)ipk(t,y) defines the kernel of an operator R(t) which 
solves the von-Neumann-Heisenberg equation with the corresponding initial 
datum, and vice versa. 

The Schrodinger-Poisson picture is more suitable for our present pur- 
poses, which are as follows: We investigate the nonlinear stability of certain 
steady states of the Schrodinger-Poisson system, i. e., of solutions of the form 
ij)k(t,x) — e l ^ kt (f)k(x) with energy levels fik G H, and we prove the existence of 
such steady states. To our knowledge, the stability problem has not yet been 
investigated. The existence of steady states has been considered by different 
methods before, cf. il6|, O, III 



Our approach is motivated by analogous results for the Vlasov-Poisson 
system which arises as the classical limit of the Schrodinger-Poisson sys- 
tem. Both systems share the following property: The total energy of the 
system is conserved along solutions — indeed, the dynamics can be inter- 
preted as the "Hamiltonian flow" induced by the energy functional — , but 
the steady states are not critical points of the energy. On the other hand, 
there exist additional conserved quantities, the so-called Casimir functionals 
[§], such that a given steady state is a critical point for the appropriately 
chosen energy-plus- Casimir functional Tic- The energy-Casimir method was 
first used to prove genuine, nonlinear stability for fluid-flow problems by 
Arnol'd in the 1960's, cf. |2], [|. Some of the background of this method 



can be found in [15]. More recently, the energy-Casimir method was adapted 



to problems in kinetic theory, in particular the Vlasov-Poisson system, cf. 



[0, [11], |T|, 0, |I§ 0, |2T], g|. When applying this method there is a 
sharp contrast between the plasma physics situation and the stellar dynam- 
ics one, where the sign in the Poisson equation is reversed: The quadratic 
part in the expansion of the energy-Casimir functional at the steady state is 
positive definite in the plasma physics case while it is indefinite in the stellar 
dynamics case. Therefore, in the former case the method applies in a straight 
forward manner, cf. [20|, while in the latter careful investigation of 

the behavior of the energy-Casimir functional along minimizing sequences 
is needed and leads to nonlinear stability only for such steady states which 
are obtained as minimizers of this functional. The present paper addresses 
the plasma-physics case for the Schrodinger-Poisson system, and thus the 
approach should be more like the former case for the Vlasov-Poisson system. 

This is indeed so: In Section 3 we show that steady states (^o^o) 
from a specified class are nonlinearly stable, and we do so by estimating 
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Hc(4>,X) — 7~Lc(i J o, Ao) from below by an expression which is quadratic in 
(■0, A) — (^o, Ao), where (if>,\) is some other, 'close-by' state. In Section 4 
we construct a functional which is in some sense dual to a given energy- 
Casimir functional. As shown in Section 5 this dual functional has a unique 
maximizer, which is a steady state, and nonlinearly stable by Section 3. We 
emphasize that — as opposed to the stellar-dynamics situation for the Vlasov- 
Poisson system — the stability analysis and the existence analysis are inde- 
pendent from each other; the connecting Section 4 is included to put both 
parts into a common perspective. Before going into all this we introduce the 
class of steady states respectively Casimir functionals under consideration, 
derive some preliminary estimates, and fix some notation. 



2 Preliminaries 

As state space for the Schrodinger-Poisson system we use the set 

S:={(ip,x) | v = (V> fc )fc€]N c H*(n) n H 2 (n) 

is a complete orthonormal system in L 2 (f2), 
A = (A fc ) fceN G I 1 with A fc > 0, fee IN, 

E fc Afc/l A^| 2 <oo}; 

J2k always means Y^kLi- O ur notation for the Sobolev spaces H 2 and Hq is 
the standard one; by || • \\ p we will denote the norm in the usual IP space. 
For (ip, A) G S we have 

and V^a denotes the Coulomb potential induced by n^ } \, i. e., 
A^,,a = -n^A on Q, V^ x = on 90; 

note that I^a G ffo flif 2 (f2) by the energy bound and Sobolev inequali- 
ties. For every initial state (ip(0),X) G S there is a unique strong solution 
[0,oo[9t^(t) of O)-(0) with (V'(*) ) A)G<S, cf. [g. Throughout the 



paper, potentials V are real-valued while quantum states ipk are complex- 
valued. The energy of a state (^,A) G<S is defined as 

W,A) :=E fc ^ / |V^r + ^/n^,A 
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£ fc A fe /|V^r + ^/|VV^ A | 2 ; 



integrals always extend over the set Q. The energy is conserved along solu- 
tions of the Schrodinger-Poisson system , indeed, the system (|1.1|)-(|L4|) can 
be written in the form 



dt 



0. 



where the bar denotes complex conjugation. 

To assess the stability of a given steady state we employ an energy- Casimir 
functional 

W c (^A):=£/7(A fc )+W,A) 

with the real-valued function C defined appropriately. Clearly, Tic is a con- 
served quantity for the Schrodinger-Poisson system. 

The class of functions which generate the Casimir functionals will now be 
specified: We say that a function / : 1R — > 1R is of Casimir class C iff it has 
the following properties: 

(i) / is continuous with f(s) > 0, s < So and f(s) = 0, s > So for some So G 
]0,oo], 

(ii) / is strictly decreasing on ] — oo,so] with lim s _ + _ 00 /(s) = oo, 

(iii) there exist constants e > and C > such that 

f{s)<C(l + s)- 7/2 - e , s>0. 

For feC, 

f(a) da, seIR, (2.1) 

defines a decreasing, continuously differentiable, and non-negative function 
which is strictly convex on its support, and 

F{s)<C{l + s)- 5/2 - e , s>0. 
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In passing we note that by adjusting various exponents our results easily 
extend to general space dimensions. 

Remark 1 (a) A typical example for f EC is the Boltzmann distribution 
/(s) = e _/3s with f3>0, where the cut-off level sq = oo. Another exam- 
ple, which also decays exponentially for s-^oo, is given by the Fermi- 
Dirac statistics: 

where C > and e > are positive parameters. 

A function / with /(s)=0 for s>s with s GlR will yield a steady 
state consisting of a finite number of quantum oscillators. 

(b) We could generalize the assumption (iii) to requiring that both /(— A + 
V) and F(— A + V) are of trace class for (smooth) potentials V>0, cf. 
Lemma [I] (b) below. However, we prefer to keep our assumptions on / 
explicit. 

Lemma 1 Let f EC. 

(a) For every (3> 1 there exists C = C{(3) E IR such that 

F(s)>-(3s + C, s<0. 

(b) Let VeHq(Q) be non-negative on Q. Then both f(—A + V) and 
F(— A + V) are trace class. 

Proof. Part (a) is straight forward from assumption (ii) and the definition of 
F. As to (b), let (fik) denote the sequence of eigenvalues of — A+V. Then, 
since V is non-negative and F decreasing, 

where /i° denote the eigenvalues of —A. For the latter we have the well- 
known estimate that the number of such eigenvalues less than some [i G IR 
grows like /z 3//2 for /x— >oo, which implies that the right hand sum is finite, 
and F(— A+V) is trace class. Since / decays faster than F the same holds 
true for f{—A+V). □ 
At several points the following technical observation will be useful: 
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Lemma 2 For i/j G H£(Q) nH 2 (n) with ||^||2 = 1 andVeH^Q), V>0, we 
have 

F((tl>, (- A +V») < (V, F(- A +V» 
with equality ifip is an eigenstate of — A+V. 

Proof. Denoting the spectral measure associated with — A +V and if) by da 
the claim translates into the inequality 



f(^J ada^j < J F(a)da 



which holds due to the convexity of F and Jensen's inequality. □ 
To conclude this section we make precise the notion of a steady state of the 
Schrodinger-Poisson system: A quadruple (ipo, Ao,//oj Vo) with (-00, Ao) GS, 
A*o — (A*o,fc) £ an d Vo G H 2 (VL) fl i?o (^) is a steady state of the Schrodinger- 
Poisson system ^l.l\ )- ftT7%) iff 

(- A+V )ip 0jk = Ho,k?Po,k, fee IN, (2.2) 

and 

A V = -n = -]T fc Ao, fc |<M 2 , (2.3) 

where the energy levels /i ,fc and occupation probabilities Ao,fc are related 
through an equation of state of the form 

Ao,fc = /W), fee IN, (2.4) 

with some / GC 



Remark 2 If (-0 O , Ao,/^oj^o) satisfies the equations Q2.2|), (|2.3|), ( |2.4| ) with 
/ G C then the estimate 

£ fc A ,fc||V'0,*||j?a<OO. 

follows and thus in particular (if;, A) G 5. To see this we use ( |2.2|) and estimate 

E fc Vllv^o, fc |l2+ / |vy r=E fc ^,fc/(^o,*)^ cr E fc ( 1 +^.*)" (5/2+e) < 00 

by assumption (iii) on / and the asymptotic behaviour of fio,k- Thus, by the 
Sobolev inequality, 

1Kb < EtVfcll^/VHe < °°> 
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and V eW 2 ' 3 (Q,)cL°°(£l) follows. Again from ( pT2|) we conclude that 

<c(i + Y, k (i+^ k )- {3/2+ * ) )<°o. 



Given / e C we still need to specify the corresponding Casimir functional: 
With F given by Q2.1|) , its Legendre or Fenchel transform is defined by 

F*(s):=sup(As-F(A)), seIR, (2.5) 
AeR 

and the energy-Casimir functional corresponding to / is 

Hc(ip,X):^J2 k F *(- x k)+n(iP,X), (ip,X)eS. (2.6) 
Note that since F' = —f has an inverse on ] — oo,so[, 

F*(s) = fj-\°)do (2.7) 

for — oo = — /(— oo) < s < 0, and all —\ lie in this interval. 

Obviously, only the values of / G C on the interval ]0, oo[ are significant for 
the following theory. However, for technical reasons we consider the functions 
/ defining the equations of state as defined on all of H. 



3 Stability 



In the present section we shall establish the following result on nonlinear 
stability: 

Theorem 1 Let (ipo, \q,/iq, Vq) be a steady state of the Schrddinger-Poisson 
system with 

A ,fc = /(A i o,fc), fee IN, 

for some fEC, and (i[) ,\ ) G«S. Then this steady state is nonlinearly sta- 
ble in the following sense: If £i— » (tp(t),X) is a solution of the Schrddinger- 
Poisson system with initial datum (^(0),A) G<S then 



i>{t),\- 



Wn 



<Wc(^(0),A)-7Yc(Vo,A ), *>0. 
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We recall that Ttc is defined by ( |2.6|) for the given function / and note that, 
clearly, the right hand side in the estimate above becomes arbitrarily small 
if (ip(0),X) is close to (if) ,\ ) in the appropriate topology. The main step in 
the proof of Theorem [I] is to show the following estimate: 

Lemma 3 LetVeH%(Sl), V>0. Then 



F*(-X k ) + X k 



lW> fc | 2 + V>* 



>-Tr[F(-A+y)], {iP,\)eS, 



with equality for (ijj, A) = (ipv, Ay), where ^y = (/0y&) e -^o(^) is an or- 
thonormal sequence of eig en functions of — A +V with eigenvalues \iy = 

(jlV,k)> and Ay=(Ay fe ) = (/(>yfc))- 

Proof. The fact that F and F* are related by conjugacy implies that 



F*(-A) + A/i> M[F*(-s)+ Sf ,} 
sen 



-sup[- F*(s) + sfj] = -F**(n) 
sen 

-F(fi), A,/xeIR. 



(3.1) 



We substitute At for A and 



Wk\ 2 + V\il> k \ 2 =<Vfc,(-A+V> fc ) 



for /i and sum over k to find 



F*(-\ k ) + \ h 



|W> fc | 2 + V> fc | 



>-J2 k F((^,(-A+V)^ k )) 

>-J2 k (^F(-A+V)^ k )) 
= — Tr[F(— A+V)] 



by Lemma and the definition of trace. 

Now suppose that (if), A) = (^y, Ay). Since by definition each ipv,k is an 
eigenfunction of — A+V the /i^ defined above are the corresponding eigen- 
values /iy*;, and 

Tr[F(-A+V)}=J2 k F M- 

On the other hand we have Ay*, = f{nv,k) — ~F'(^ Vjk ) which by conjugacy is 
equivalent to fiv,k = F*'(— Ay&), fee IN. This implies that 



J2 k F (^v,k) = -J2 k l F *(- x v,k) + \v,kw,k 



and the proof is complete. 



□ 
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Remark 3 In Lemma |3| equality holds if and only if (if>,\) = (ipvAv)- This 
follows from the strict convexity of F, but we make no use of this observation 
in the rest of the paper. 

Proof of Theorem |Tj. Let V = V^ a be the potential induced by (if), A) £ S. 
Then 

kw-Woll* 



~/iwr+/A^o+i/iw i 



= Hc{if>,\)- 
= Hc{if>,\)- 



Y, k ( F *(-h)+h J |v^r)-^/|vy | 2 - / 



V F*(-A fc ) + A fc / [|W> fc | 2 +Wfc 



afv 

IVV0I 2 



-Tr[F(-A+F )]-i/ 1 Wo 



£ fc (V(-A ,,) + A , fe J (|W> ,fc| 2 + Vbl^o^l 2 )) - -/ 



2 / |W °' 



■ Hc(i),\)-H c (i)oAo) 



where we have used Lemma |3] twice. Given a solution with (ip(0), A) G S 
we may substitute (if>(t),\) G5 into this estimate, and since He is constant 
along solutions the assertion follows. □ 



4 Dual functionals 



Our aim for the rest of this paper is to prove the existence of steady states 
which satisfy the assumption of our stability result. For each / eC a corre- 
sponding steady state will be obtained as the unique maximizer of an ap- 
propriately defined functional. In the present section we derive this dual 
functional from the energy- Casimir functional used in the stability analysis. 
The relation between these functionals is of interest in itself, but it is not 
used in the proofs of our results. Throughout this section we fix an element 
/ G C. We move to the dual functional in two steps. First we apply the saddle 
point principle and define, for A > fixed, 

1 



£(V,A,V»:=V 



F*(-\ k ) + \ k / |V^| 2 + m 



- w 
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where ip=(ipk) is again an orthonormal system in L 2 (Q), \El]_ = {(cr fe ) G 
^Icfe^O, fee IN}, and KGff^fi) may now vary independently of ijj and A. 
The role of the parameter o E IR (Lagrange multiplier) will become clear 
shortly; the relation between Tic an d this new functional is as follows: 

Remark 4 For any tp,X,a, 



and the supremum is attained at V = V$ t \. In fact, integration by parts and 
some computations show that 



As second step on our way to a dual variational formulation we reduce the 
functional Q to a functional of V and a as follows: 



where the infimum is taken over all A G l\ and all orthonormal sequences i/j 
in L 2 {Vt). It is this functional which will have a unique maximizer in the next 
section, which is then a steady state. First however, we need to bring it into 
a different form: 

Remark 5 The infimum in the definition of $ is attained at ip = {ipv,k)i an 
orthonormal sequence of eigenstates of — A+V with corresponding eigenval- 
ues and A = Ay where Xy,k = f(^v± + &), fee IN. Moreover, 



sup£(V,A,V>)=H c (^,A)+a[£ fc A fc -A , 



(4.1) 



g^,X,V,a)=H c ^,X)+ ( r[£ k \ k -A\-l\\VV^-VV 



$(V»:=inf0ftM,V» 



(4.2) 




To see this, recall Lemma [3] and Remark [3] and observe that f {• -\- a) EC for 
any a G IR, provided f EC. 
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5 Existence of steady states 



In the present section we shall for each state relation f EC and each total 
charge A>0 construct a unique maximizer of the functional $, which is 
then a steady state of the Schrodinger-Poisson system. We consider only 
non-negative potentials and use the notation 

Hl + (p):={veH*(n)\v>o}. 

Theorem 2 Let f EC and A > be given. The functional 

<5>:Hl + (Q)xlR3(V,o-)^-^J |W| 2 -Tr [F(-A + V + a)\ -a A 

is continuous, strictly concave, bounded from above, and coercive. In partic- 
ular, there exists a unique maximizer (Vo,o~o) o/$. If we define ipo — (ipo,k) 
as the orthonormal sequence of eigenstates of the operator — A+Vq with 
corresponding eigenvalues fio,k and Ao,fc := f(^o,k + o~o), then (ip , Ao,jUo, Vo) 
is a steady state of the Schrodinger-Poisson system with J2k^o,k = A and 
(ipo,Xo)ES. 

Note that o~q plays the role of a (constant) Fermi level here. 
Proof. $ is strictly concave: The first term of <3> is evidently concave. To 
show the strict concavity of the second term i.e., the strict convexity of 
Tr[F(-A + V + a)], let (Vj,o-j) E Hq + x H, j = 1,2, ae]0,l[, and (j)EH 2 n 
Hq. By convexity of F and Lemma 0, 

F(( ( f),a(-A+V 1 + a 1 )(f ) +(l-a)(-A+V2 + o- 2 )(f ) }) 

<a((l>,F{-A+V 1 + a 1 )<f>) + (l-a)((l>,F{-A+V 2 + <T 2 )<f>). 

Now we substitute ipk for (ft, (ipk) an orthonormal sequence of eigenstates 
of a(— A +V\ +<Ti) + (1 — Oi)(— A +V2 + cr 2 ), and sum over k to obtain the 
convexity estimate for Tr [F(—A + V + a] . If we have equality in this estimate 
then 

(ip k ,F(- A +Vi + rrO^fc) = (ip k ,F(- A +V 2 + <r 2 )if> k ), k E IN 
and thus V\ — V 2 and <j\ — <t 2 . 
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$ is bounded from above and coercive: Since F is non-negative, the critical 
case in the coercivity estimate is a < 0. Let fi y denote the ground state energy 
of — A + V with corresponding ground state if) . Since F is non-negative and 
satisfies estimate (a) in Lemma p] we have for a < —fly, 

HV,a) < ~J\VV\ 2 -{t v ,F(-A+V+a)iPy)-aA 
= _l||VF| 2 -F(-/, v , + a)-aA 
<-lf\VV\ 2 + (P-A)a + P l x v -C, 



where we choose (3 > A. Also 



Hv = , eH M, , / -!^| 2 + V|0| 2 : —f Y<C\\\Y\\ U ,, 



1 

volfi 



choosing (f>: = 1/vvolO. Together with the estimate above and Poincare's 
inequality this implies that for a < —Ci\\V\\ H i we have 



^(V > a)<-C 2 \\V\\ 2 Hh +C 3 \\V\\ H x + {p-A)a + C A (5.1) 



where the constants C\,Ci,C^,Gi are positive and f3>A, cf. Lemma |1| (a). 
On the other hand, by the non-negativity of F and Poincare's inequality, 



§{V,a)<-C 2 \\V\\ 2 H i-aA, (5.2) 



and (|5.1| ) and ( p.2|) together imply that $ is bounded from above and coercive. 
Existence of a unique maximizer: The existence of a unique maximizer of 
$ is standard, cf. for example || Ch. II, Prop. 1.2], provided $ is upper 
semi-continuous. This in turn follows from the fact that $ is concave and 
bounded from below, at least locally, cf. @, Ch. I, Lemma 2.1]: The only 
term for which this may not be immediately obvious is the trace term, but 

Tr[F(-A+V + a)]<Y; k F(Hk + vo)<oo 

where fif. are the eigenvalues of — A +V and a > <t for arbitrary o~ 6 IR. 
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(ipo, \),A*o> Vo) i s a steady state: Since F' = —f, the stationarity of $(Vo,cr) 
with respect to a implies 



■ Tr[f(-A+V + a ))-A 



da 

= J2 k f(^o,k + (T ) - A = ^2 h Xo,k - A 

so that I]fc-Vfc = A as claimed. In order that (ipo, Ao,A*o,^o) is a steady state 
we need to show that 

A ^o + E fe Vcl<M 2 = 0. (5.3) 

To verify this we observe that Vq, being a maximizer of $(-,o"o), satisfies the 
Euler-Lagrange equation 

AV (x)+K f{ _ A+Vo+(7o) (x,x) = 0, xett, (5.4) 

where Ki is the kernel associated with a trace-class operator L. In our case 

Kf(- A+ V +a )(x } x)=Y, k f(^0,k + ^o)\'ipO,k\ 2 ( X ) ( 5 - 5 ) 

and ( |5.3| ) follows from (|5T^), (|5.5|), and the fact that by definition, X 0jk = 
f(^o,k + o"o). As to the proof for (tp ,\ ) 6 5 we refer to Remark ||. □ 
In view of the relations between our various functionals derived in the pre- 
vious section it is of interest to note: 

Remark 6 If (Vo,<t ) is the maximizer obtained in Theorem |2| and 
(■00, Ao,//q,Vq) is the corresponding steady state, then 

$(y 0) (7o)=W C (V'0 ) Ao). 

To see this, note that by ( |4.2[ ) we have 

<5>{V ,a ) = g{tp ,\ ,V ,o )<H c {4JoAo), 

where equality holds iff V is the maximizer of Q(ip ,X ,V,a ) on Hq] note 
that here Q is independent of a since Y^k -Vfc = 1- This, on the other hand, is 
equivalent to the fact that Vq is the solution of the Poisson equation ( |5.3[) . 



14 



References 



[1] Angelescu, N., Pulvirenti, M., Teta, A.: Derivation and clas- 
sical limit of the mean-field equation for a quantum Coulomb system: 
Maxwell-Boltzmann statistics. J. Stat. Phys., 74, 147-165 (1994) 

[2] Arnol'd, V. I.: Conditions for nonlinear stability of stationary plane 
curvilinear flows of an ideal fluid. Sov. Math. Dokl. 6, 773-776 (1965) 

[3] Arnol'd, V. I.: On an a priori estimate in the theory of hydrodynam- 
ical stability. Amer. Math. Soc. Transl, Series 2, 79, 267-269 (1969) 

[4] Brezzi, F., Markowich, P. A.: The three-dimensional Wigner- 
Poisson problem: existence, uniqueness and approximation. Math. Meth. 
Appl. Set. 14, 35-62 (1991) 

[5] Casimir, H. G. B.: Uber die Konstruktion einer zu den irreduziblen 
Darstellungen halbeinfacher kontinuierlicher Gruppen gehorigen Differ- 
entialgleichung. Proc. R. Soc. Amsterdam, 34, 844-846 (1931) 

[6] CASTELLA, F.: L 2 -solutions to the Schrodinger-Poisson system: exis- 
tence, uniqueness, time behaviour, and smoothing effects. Math. Mod. 
Meth. Appl. Set. 7, 1051-1083 (1997) 

[7] DOLBEAULT, J., MARKOWICH, P. A., Unterreiter, A: On singular 
limits of mean-field equations. Arch. Rat. Mech. Anal, to appear 

[8] Ekeland, I., Temam, R.: Convex Analysis and Variational Problems. 
New York: North-Holland Publishing Company, 1976 

[9] Gasser, I., Illner, R., Markowich, P. A., & Schmeiser, C.: 
Semiclassical, t^oo asymptotics and dispersive effects for the Hartree- 
Fock systems. Math. Mod. and Num. Anal. 32, 699-713 (1998) 

[10] Guo, Y.: Variational method in polytropic galaxies. Arch. Rational 
Mech. Anal, 150, 209-224 (1999) 

[11] Guo, Y.: On the generalized Antonov's stability criterion. Contem. 
Math. 263, 85-107 (2000) 



15 



[12] Guo, Y., Rein, G.: Stable steady states in stellar dynamics. Arch. 
Rational Mech. Anal. 147, 225-243 (1999) 

[13] Guo, Y., Rein, G.: Existence and stability of Camm type steady 
states in galactic dynamics. Indiana University Math. J., 48, 1237-1255 
(1999) 

[14] Guo, Y., Rein, G.: Isotropic steady states in galactic dynamics. Com- 
mun. Math. Phys., to appear 

[15] Holm, D. D., Marsden, J. E., Ratiu, T., & Weinstein, A.: 
Nonlinear stability of fluid and plasma equilibria. Physics Reports, 123, 
Nos. 1 and 2, 1-116 (1985) 

[16] Markowich, P. A.: Boltzmann distributed quantum steady states 
and their classical limit. Forum Math., 6, 1-33 (1994) 

[17] NlER, F.: A stationary Schrodinger-Poisson system arising from the 
modelling of electronic devices. Forum Math. 2, 489-510 (1990) 

[18] NlER, F.: A variational formulation of Schrodinger-Poisson systems in 
dimension d< 3. Commun. Partial Differential Equations 18, 1125-1147 
(1993) 

[19] NlER, F.: Schrodinger-Poisson systems in dimension <i<3: the whole- 
space case. Proc. Roy. Soc. Edinburgh Sect. A 123, 1179-1201 (1993) 

[20] Rein, G.: Nonlinear stability for the Vlasov-Poisson system - - the 
energy-Casimir method. Math. Meth. Appl. Sci, 17 , 1129-1140 (1994) 

[21] Rein, G.: Flat steady states in stellar dynamics — existence and stabil- 
ity. Commun. Math. Phys. 205, 229-247 (1999) 

[22] Wolansky, G.: On nonlinear stability of polytropic galaxies. Ann. 
Inst. Henri Poincare, 16, 15-48 (1999) 

[23] Wolansky, G.: Static solutions of the Vlasov- Einstein system. 
Preprint, 2000 



16 



